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SUMMARY 


The vibration eigenvalues and modes have been analyti- 
cally determined for the helicoidal shell and its limiting 
case, the skewed flat plate, A general analysis is given 
wherein the energy equation employed is derived from 
Hamilton's Law of Varying Action, No numerical data has 
been found in the literature for the helicoidal shell. Re- 
sults for spherical caps and complete conical shells, ob- 
tained by the method of this paper, are compared to exact 
solutions from the literature for these two shell geometries. 
The authors consider this paper to be a significant step 
toward an analytical solution to the thermally stressed, 
variable thickness, turbine blade immersed in both a fluid 
flow field and in an acceleration field. 
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INTRODUCTION 


In recent years emphasis in solid mechanics has shifted 
from plates to shells. This can be attributed to the fact 
that thin plate vibrations have been widely explored for 
simple geometries such as the circle, the rectangle, and the 
rhombus [1], Exact solutions were found for isolated cases. 
However, for many boundary conditions the analyst was (and 
is) forced to seek an approximate solution. "Approximate" 
usually meant a Rayleigh-Ritz type solution. "Beam" func- 
tions became popular among plate analysts who attempted 
solutions by Rayleigh-Ritz since they already satisfied the 
geometric boundary conditions. The results in many cases 
were less than satisfactory. This inaccuracy was generally 
attributed to the fact that the method was "approximate". 

With the advent of finite element analysis, the analyst 
had at his disposal an approximate method which took much 
less intuition than Rayleigh-Ritz (the Rayleigh-Ritz method 
requires one to assume a solution in the form of "admissible 
functions). A predominance of the finite element technique 
has come about for numerous reasons, but two major factors 
are the easy satisfaction of the boundary conditions and the 
ease of handling complicated geometric shapes. However, us- 
ing concepts long understood for satisfaction of boundary 
conditions [2], the Ritz method can now be applied to geo- 
metrically complicated structures with great accuracy for 
both stationary and non- stationary problems. 

Since finite element analysis is now widely accepted as 
the best way to solve the difficult problems, the Rayleigh- 
Ritz procedure has been pushed aside [3]. In the authors' 
opinion, it is unfortunate that analysts always think of 
Rayleigh-Ritz instead of just Ritz [4]. Rayleigh's intent 
was that of choosing one function, which if at all close to 
the true mode shape, would yield a fair, if not good, 
approximation for the frequency by means of Rayleigh's quo- 
tient. However, Ritz's idea was that of choosing a series 
of suitable approximating functions. Ritz [4] showed that 
such a series, with a finite number of terms, converges to 
an upper bound of the true solution. His proof, however, and 
the resulting direct solutions apply only to stationary 
systems. Using concepts [5] associated with Hamilton's Law 
of Varying Action [6] the authors have obtained direct solu- 
tions to non-stationary and/or non-conservative systems. 

When calculating the transient motion of a deformable body, 
a major simplification for achieving direct solutions (as 
with the equations of Lagrange) results when one uses the 
free vibration modes and normal coordinates, whether or not 
this procedure uncouples the system. Thus, the generation 
of free vibration modes and frequencies continues to be of 
great importance in the solution of the problems of the 


displacement or motion of physical systems. As a result, 
the authors have restricted this paper to a particularly 
difficult (and unsolved to the author's knowledge) but con- 
servative, stationary problem: the free vibration modes and 
frequencies of a helicoidal shell. 

A fact which seems to have been overlooked, possibly 
because of the concepts associated with the Rayleigh-Ritz 
method, is the accuracy obtainable when power series are 
taken as the "admissible" functions. In this paper. Appendix 
I, and in other work, the degree of accuracy which can be ob- 
tained by direct solutions of the energy equation has been 
compared by the authors to exact solutions where available. 

In all applications (particles, rigid bodies, beams, plates, 
and shells), when finite time intervals and finite space 
intervals are properly employed, the degree of accuracy far 
surpassed expectations. The Ritz method, in every case, con- 
verged to the exact solution, for the constraints assumed, 
in the same sense that any number representing the solution 
to a transcendental equation can be considered to be exact. 
The results of these endeavors all indicate that the idea of 
"approximation" associated with the Rayleigh concept should 
not be associated with the Ritz method. 

As mentioned above, this study is restricted to the 
stationary problem of simple harmonic vibrations. For the 
helicoidal shell, the displacement field is completely des- 
cribed by the displacements u, v, and w, all three of which 
are functions of the orthogonal surface coordinates x and y. 
In the limiting case of the flat plate, only w is considered, 
where w is a function of the rectangular inplane coordinates 
x and y. Inplane (or membrane) vibration for the flat plate 
will not be treated, although the authors have completed 
such solutions for the rectangular cantilever plate [2]. As 
has been shown for a completely free elliptical plate [7]> 
and for plates with various other boundary conditions (Refs. 
[2], [8], when power series are used the Ritz method yields 
very accurate mode shapes as well as accurate frequencies. 


This work was supported in part by the Department of 
Aeronautical and Astronautical Engineering, The Ohio State 
University, and in part by the NASA Langley Research Center 
through Grants NGL-36-008-109 and NOR-36-008-197. 
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CHAPTER I. 


THEORY 

Hamilton ’ s Law 


In attempting an energy solution, it seems that the 
analyst is always confronted with which energy equation or 
theorem to use as a starting point. One form which embodies 
all other energy equations of the mechanics of solids, is 
that enunciated in 1834 by Sir William Rowan Hamilton [6] 
and called by him the "Law of Varying Action". It is not 
the well known Hamilton's Principle. Indeed Hamilton's 
Principle is a special case of Hamilton's Law (throughout 
the text Hamilton's Law will refer to the Law of Varying 
Action). Hamilton's Principle is usually written [9] as, 

<$ f (T + W) dt = 0 (1) 


where T is the total kinetic energy of the system and W is 
the total work of forces acting on or within the system. 
Hamilton's Law, which is much more general [6], is expressed 
as. 


6 / (T + W) dt - 
t n 9q i 


^1 

6q ± | = 0 

t 0 


( 2 ) 


In the variational calculus, a stationary solution is, by 
definition, a solution of a given system for which 




0. Eq. (2), which reduces to Eq. (1) for 


stationary motion, applies to non-stationary motion as well 
[5]. No matter what the motion is, Eq. (2) will yield either 
a direct solution (as in this paper) or the differential 
equations of the system. 


Strain-Displacement Equations 

It is assumed that the material of the shell is isotropic 
and obeys a linear stress-strain relationship. In addition, 
it is assumed that the oscillations cause infinitesimal dis- 
placements which can be described by a set of linear strain- 
displacement relations. These strain-displacement relations 
come about as a result of considering the deformation of the 
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middle surface. The additional assumption that normals to 
the middle surface remain normal throughout the deformation 
process, always made in thin shell theory, then locates all 
other points in the body relative to the middle surface. 
Since differential geometry is not the topic of this paper 
the authors will merely restate the linear strain- 

displacement relations for thin shells given by Novozhilov 

[10]. 


e 


x 


1 3U , 

A 3x 


J_ M v 
AB 3y 


+ 



y B 3y AB 3x R y 


e = ? ~ (v/B) + A J_ (u/A) 
xy A 3x B 3y 


A 3 , 


K. 


_-l _3_ /I 3w _ _u_v _ _1_ _3A .1 
A 3X A 3X 13 ' ^ 


FL 


3w 

AB 3y V B 3y 


£ 


K = - i JL ri — - — ^ 1 3B .1 3w _ U ■ 

y B 3y m 3y R. ' “ AB 9x '■A 3x R x ' 

u 


V - 2 ( 3 2 w 1 3A 3w 1 3B 3w 

*xy AB -3x3y A 3y 3x B 3x 3y' 


+ JL rl M _ JL M T7% 4 . JLri iL _ JL 

R„ ^B 3y AB 3y u; + R t /A 3 x AB 3 x v ^ (3) 

-*• y 


For any point on the surface of a right circular cylinder, 
where x is the axis of the cylinder and y is the circumfer- 
ential angle measured from the vertical (clockwise positive, 
see Fig. 3) 5 the following geometrical relations are true: 

"R =2 OO 

n x 

R y = R = constant 

A = 1 

B = R = constant 
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8A _ d_A _ JIB _ .3B _ \ _ 3 r-o \ JL/ R 'i = JL/'p \ 

3y 3x 3y 3x 3y^ x ; “ Sx^x' Sy^^V' 3x^ K y' = 0 (4) 


where R is the radius of the cylinder, 
cations Eq. (3) becomes. 


With these simplifi- 


e = M 
x 3x 


e = I (1Z + w) 
y R v 3y ; 


_ 3v , 1 3u 
'xy ax R 3y 


K x = 


3 2 W 

3 x 2 


K y = 


K 


xy 


2- - 

1 3 w _1_ 3v; 

r 2 3 y 2 R 2 3y 

2 - - 

- 2 3 w 2 3v 

R 3x3y R 3x 


(5) 


The total strain of any point in the body is given by. 


e x = e x + z K 5 


e y e y + z K y 


e xy E xy + z K xy 


( 6 ) 


where z is the distance from the middle surface. For the 
plane stress problem (that is a z is negHigible, thin shell 
assumption) the stress at any point in the body is given by. 


c x (i— (e x + v e y) 


°y “ (• ^ e y + v e x^ 

g — 

CT xy = 2 ( l+v ) e xy 


(7) 
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Work 


For an elastic body, the work, W, is that due to the 
internal elastic restoring forces, the internal or "body" 
forces, and the external or surface forces: 

W = - fffCa x e x + o y e y + a xy e xy ) R dx dy dz 

+ ///( u + Fgy v + Fg z w) R dx dy dz 

+ //(F g x u + F Sy v + Fg z w) dS (8) 

Taking the variation of Eq. (8) with respect to the displace 
ments (as in the principle of virtual work) yields, 

6W = - ///(o x 6e x + ay 6e y + a xy <$e X y) R dx dy dz 

+ fff ( Fg 6u + Fg 6v + F B 6w) R dx dy dz 

x y 2 

+ //(Fg 6u + Fg 6v + Fg 6w) dS (9) 

x y 2 


Kinetic Energy 


Rotary inertia is negligible for thin plate and thin shell 
vibrations, thus, neglecting rotary inertia, the total 
kinetic energy is, 

T = Tj- /// p (u 2 + v 2 + w 2 ) R dx dy dz (10) 

and 



Iff p(u 6u + v Sv + w 6w) R dxdydz | (11) 


Taking the variation of Eq. (10) and substracting Eq. (11) 
in Hamilton’s Law, 
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r 


t n 5 m 2 _ 

{ <T at - 84T eq i i = 
t 0 t 0 


r- ^ _•_ -i 

/ ///p(u 6u + v 6v + w Sw) R dx dy dz dt 

to L 


. . « 1 1 
- ///p(u 6u + v <5v + w 6w)R dx dy dz | (12) 

t 0 


Integrating the first term on the right hand side of Eq. (12) 
by parts with respect to time yields. 


/ 6T dt 

to 


f- *<.! I 1 

S Q1 t 0 


-h 


fff p(u 6u + v 6v + w 6w)R dx dy dz 


dt 


(13) 


Energy Equation 

Substituting Eqs. (9) and (13) into Eq. (12), the energy 
equation (Hamilton's Law) becomes, 

tn r ui. si. 

/ -/// p(u 6u + v 6v + w <Sw)R dx dy dz 

to L 

-fff(o x 6e x + Oy 6ey + F X y 6e X y)R dx dy dz 

+///(F R 6u + F r Sv + P B Sw) R dx dy dz 
B x °y D z 

+ //(P Sx <5u + F Sy Sv + F Sz 6w)dS J dt = 0 (14) 
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Eq. (14) applies to non-stationary as well as stationary 
motion. For the case of no external' forces acting, assume 
simple harmonic motion (sin tut or cos tut or e 10 ^ will all 
yield exactly the same result in this stationary problem) 

u(x,y,t) = u (x,y ) cos tut 

v(x,y , t ) = v(x,y) cos tut 

w(x,y,t) = w(x,y) cos cut (15) 

Since the strain-displacement relations are linear and the 
body is elastic (thus ^(xjyjt) = a x (x,y)cos tut, e x (x,y,t) = 
e x (x,y)cos tot, and likewise for cT y , ay y , e y , e xy ), Eq. (14) 
becomes, upon substitution of Eq. (15), 


|^tu 2 ///p(u 6u + v <5v + w 6w)R dx dy dz 
- fff(a x 6e x + a y 6e y + a xy 6e X y)R dx dy dz 
t l 2 

f cos tut dt = 0 


(16) 


Since the integral on time in Eq. (16) does not vanish for 
any tq greater than to, one must have, 

2 

tu ///p(u 6u + v 6v + w 6w)R dx dy dz 


— f S f ( p x 6e x t cty 6e y t o X y 6e xy )R dx dy dz — 0 (17) 


Eq. (17) is the familiar equation of the principle of minimum 
potential energy which is applicable only to stationary systems. 
If the structure is symmetric about the middle surface, the 
neutral surface and middle surface are coincident. Assuming 
such is the case, upon substitution of Eqs. (6) and (7) 3 Eq. 

(17) becomes. 
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to 



yi 

/ 


x o yo 


ph(u <Su + v Sv + w Sw)R dy dx 


(1-v 2 ) 


x i yi 
/ / [h 

x o yo 


{(■ 


v e. 


■ )6e x + ( e v + v e x )6e. 


+ Or 1 e xy «>xyl + ll K k x + » Ky)«K x 


+ (Ky + V K x )6Ky + K xy 6K X y} ]R dy dx = 0 (18) 

Let C = Eh/(l-v 2 ) and D = Eh3/12 (1-v 2 ) . Substitution of 
Eq. (5) into Eq. (18) yields three coupled equations in u, v, 
and w. 


/ f c [(u x 6u x + il=yl 
xo yo 2 


u y 6u y> 


+ (-^2^- | v x 6Uy + ^ Vy 6U X ) + (g W 6u x )]R dy dx 

o X 1 y l 

to 2 / / p h(u i$u)R dy dx = 0 

xq yo (19) 


and, 

f 1 / 1 [i ( u y < 5 v x + v u x 6 v y ) 

xo yo 

+ °(gz- v y Sv y + — 2 ~~ v x 6v x) 

+ v y 5v y + 2(l-v)v x 6v x ) 

+ 6v y) - w yy 6v y + v w xx «v y 

+ 2(l-v)w X y 6v x )]R dy dx 


— to 


2 


x i yi 

/ / ph(v 6v)R dy dx = 0 

x 0 x 0 


( 20 ) 
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and. 


f f C§(v u x 6w) + ^-(vy 6w) - v y 

x o yo 

+ V Vy Sw xx + 2(l-v)v x 6W X y) + ^-(w 5w ) 

+ D {w xx 6w xx + ^ Wyy SWyy + j^( W XX 6w yy 

+ Wyy 6w xx ) + ~^ 2 ~^ w xy <$w x y}]R dy dx 


6Wyy 


- w 


2 


/ f ph(w 6w)R dy dx = 

x o yo 


0 


( 21 ) 


Coordinate Transformation 


If the cylinder in Pig. 3 is cut along y = ir and develo- 
ped in the plane (Pig. 4), the relations between x and y and 
the skewed coordinates can readily be seen. They are. 


x 

cos a 


Rii = R y - x tan a (22) 

Prom the calculus , 

3-3 3 1 4 . 3 3~n 

dx 3x 3r T 9* 


3 _ 3 H + 3 9n 

9y 9f 9y 9^ 9y (23) 

Substituting Eq. (22) into Eq. (23), 
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3 

3x 


tan d 3 


1 _3 

COS a 3g“ R Sff 


3 _ 3 

9y 3rT 



3 2 _ 1 3 2 + tan 2 a 3 2 _ 2 tan a 3 2 

3x z COS z a 3f z R z 3^ 2 R cos a 31'3n' 

3 2 _ 1 3 2 _ tang 3 2 

3x3y " COS a 3f3r|’ R 3rf Z 

dx dy = cos a dT dr7 (24) 


Make the following non-dimensionalizations , 

t = n/y 

5 = ?/* (25) 

where y is the half-angle at £• = 0, and 1 is the slant length 
(see Pig. 4). Substitution of Eq. (25) into Eq. (24) yields, 

3 _ 1 _3 tan a 3 

3x t, cos a 3 £ yR 3n 


3 _ 1 3 

3y y 3n 

3 2 - 1 3 2 
3y 2 y 2 3n 2 


3 2 _ 1 3 2 , tan 2 a 3 2 _ 2 tan a 3 2 

3x z S, 2 cos 2 a 3£ 2 R 2 y 2 3 n 2 Ry«,cosa 3£3n 


3 2 _ 1 1_ 3 2 _ tan a 3 2 

3x3y Jlcos a y 3£3n R y 2 3 n 2 

dx dy = 5-ycosa d£ dn (26) 
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Case of the Parallelogram 


When the projection of the open circular shell in Pig. 

4 is a parallelogram (a helicoidal shell), the foregoing 
coordinate transformation will yield constant limits on the 
integrals in the energy equations [Eqs. (19) - (21)] and thus 
facilitate exact integration. Assuming a helicoidal shell 
with constant thickness h, the following equations result 
from substituting Eq. (26) into Eqs. (19)-(21), 


f f { [Ur 5Uj. + ( s in 2 a + - V - -cos 2 a) ( £/Ry) 2 u n 6u n 
1 0 [ 

- sin a (A/RyHu^ 6u n )] + (£/Ry)' 

[ cos a y ^ 6u n + v cos a 6u^ - (£/Ry) # 

(sin a) (cos a) ( v n Su^] + [(£/R) v cosa w 6u^ 

- ( £/R) ( £/Ry) v cos a sin a w Su^] 

- a> 2 g (1-v 2 ) £ 2 COS 2 a [u 6u] | d£ dn = 0 (27) 
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oa a \ijj + v eoBa u^. 8v^ 


// {(//*)pj 5*U 

- (//R-yHslna )(cosa)Il^I u,, 

+ [l + (h/R) 2 /12] [(/cosa/R r ) 2 v v 8v^] 

+ [l + (h/R) 2 /3] [(/ sina/R-y) 2 v, Cv, 

+ 6^. - (/sin a /By)(v TI C^. + v^. )J 

+ [(/cos a /R) 2 Iw 8vJ - [ ( h/R ) 2 /12y] ' 

[(//Ry) 2 (1 + (l-v)sin 2 a)w^ 6v^ + v v»^ 

- 2(//R Y ) sin a Bv^ + 2(l-v) 5^. 

- 2(l-v)(//Ry) sin a 6v^] 

- 6c) 2 | ( 1- v 2 )/ 2 coS 2 a [v 6v”|^ d£ dJ? - 0 


i 



1 1 


// £(/A) £v cos a 8w - (J?/Ry) v cosaslna 8wJ 
+ itfcosa A) 2 [ V|? 6w] - [(hA) 2 /12y] * 

[(iAy) 2 (l + (l-v)Bin 2 a) 8"^ + v v n 8w^ 

- 2 (/Ay) sina 8w^ + 2(l-v) 5w^ 

- 2(l-v)(^Ay) sina 6w^J + £(/cosa A) 2 w 6wJ 

+ jj h A ) 2 /12j jj R/i c os a ) 2 w^. 8*^ + (/Ay 2 cos a ) 2 • 
w »W 8w w + ~ (tan 2 a + vKw,^ 8w^ + ew^ ) 

+ (2(l-v) + 4 tan 2 a) Sw^ - 2(iAy) ^ * 

cos a 5w £»? + 5 w j?tj) “ 2 cos a 

^cttj COS ft 8wJ 

- 0 (29) 
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For the limiting case of the flat plate in bending, Eq. (29) 
reduces to the following; 


1 1 
/ / 
-1 0 


w 5c 6w c? + U/b) 4 w nn Sw nn 


+ ( sin 2 a + v cos 2 a.) ( «./b) 2 (w r)n 6Wg> + 6w nr) ) 

+ [2(l-v)cos 2 a + 4 sin 2 a) ()l/b) 2 (w^ r) 6Wg n ) 

- 2(*/b) 3 sina (w nn 6w ?n + w ?n 6w nn ) 

- 2 ( £/b) sina (w^g 6w ?n + w^ n Sw^) 

- (u 2 12(l-v 2 ) cos 4 a w 6w J d£ dn = 0 (30) 


where sl is the slant length of the plate and b is the half- 
width (see Fig. 2). 
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CHAPTER II 


SOLUTIONS 


Boundary Conditions 


As is well known [2], one need only satisfy geometric 
boundary conditions when using the Ritz method with dis- 
placement energy equations. For a cantilever boundary, the 
geometric boundary conditions are. 


u = 

0 

at 

x = 

0 

for 

all 

y 


V = 

0 

at 

x = 

0 

for 

all 

y 


w = 

0 

at 

x = 

0 

for 

all 

y 


w x = 

0 

at 

x = 

0 

for 

all 

y 

(3D 


From Eqs. (22) and (25) it follows that x = 0 implies 5=0. 
Eq. (26) and the last of Eq. (31) gives. 



But since w = 0 at x = 0 for all y, Wy = 0 at x = 0 for all 
y as a result. Thus in terms of the skewed coordinates 5 
and n, Eqs. (31) and (32) become, 

u = 0 at 5 = 0 

v = 0 at 5 = 0 

w = 0 at 5 = 0 

w ? = 0 at 5 = 0 (33) 
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Assumed Functions 


Using the concepts set forth in Ref. [15], direct solu- 
tions to conservative, non-conservative, stationary, and 
non-stationary systems of particles, rigid bodies, beams, 
plates, and shells have been obtained with nothing more than 
power series in Eq. (2). It has been the author's experience 
that they also give more accurate answers than beam functions 
for plate and shell vibration problems (beam functions place 
unnecessary constraints on higher derivatives at plate and 
shell boundaries for some boundary conditions). 

The geometric boundary conditions [Eq. (33)] dictate the 
following double power series. 

N M 1+ i j 

u = z Z A,- ,• 

i=0 j = 0 J 


v = Z Z B.t A ? 1+± n J 
1=0 J=0 


N M . 

w = Z Z C44 K 1+2 n J (34) 

i=0 j=0 


Eigenvalue Equations 


The 6 operator 
definitions [5], 


<5u = 


6v = 


- 3u 


3A kSL 
3 V 


9B 


k& 


6w = 


- 3w 


3C k ^ 


in Eqs. (27)- (30) now has the following 


6A kz = C k+1 


SA k i 


= S k+1 

n £ 

SB k 


6C k , - ? k+2 

n* 

6C k 

(35) 


where 6A kJl , 6B k £, and <SC k)l are arbitrary quantities (or 
"variations"), which may „be factored out of Eqs. (27)- (30). 
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The Hell co Ida! Shell 


Upon substitution of Eqs , (3*0 and (35) into Eqs, (27)- 
(29), the following coupled algebraic equations result; 


* Q { A ij ^ K1 ijkJl ~ x s + B ij K2 ijk«. 

+ c ij K 3ijk«. } = 0 


for k = 0,1. . . ,N 

% = 0,1,..,M 


MM 

± l Q K2 k*iJ + B iJ - *s M2 1JkA ] 

+ C ij K5 ijk )i } = 0 


for k = 0,1. . . ,N 
i = 0, 1. . . ,M 


(36) 


(37) 


N M , 

jf 0 { Al <5 K 3k£ij + B iJ K ^k^ij + c ij ^ K ^ijk)l 

- M3ij k ^] | = 0 for k = 0,1., . ,N 

z = 0 , 1. , , ,M (38) 


where Kl 1 j kJl , K2j_j k ^, K3ij k *, K4j_j kJl , K5j.j k £, K6ij ka , Mlij k& , 
M2ij k £ and M3ji k x, are given in Appendix II in terms of re- 
currsion formulas which result from exact integration. Eqs.. 
(36)- (38) form a set of 3 x ( N+l) x (M+l) algebraic equations in 
the (3x (N+l )x (M+l)) + 1 unknown Aji’s, By's, Cji's, and 

*s (= w (l-v 2 )£ 2 cos 2 a). Ihey are conveniently written in the following 

matrix form; 



"K1 K2 K3* 


o 

o 

s! 

) 

’a' 

< 

K2 t K4 K5 

- 4 

0 M2 0 


B 


_K3 T K5 T K6J 


o 

o 


-C. 


For a non-trivial solution, the determinant of the coefficient matrix must 
vanish. 
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The Skewed Plate 


The limiting case of the flat plate in bending, using 
Eqs. (3*0, (35) > and (30), yields the following eigenvalue 
equation: 

N M 

• ^q = 0 

for K = 0,1,2, ... ,N 

l = 0,1,2, ,M 

or in matrix form, 

{[KB] - x|[M33 } \C\ = 0 

KBijkfc is given in Appendix II. Eq. (41) is a set of (N+l) X 
(M+l) algebraic equations in the (N+l) X (M+l) + 1 unknown 
Cij ' s and X p 2 ( = io 2 ( ph/D) «. 4 cos 4 a) . Again, to have a non-trivial 
solution, the determinant of the coefficient matrix in Eq. (4l) 
must vanish. 


(40) 

(41) 
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CHAPTER III. . 


RESULTS 

The Skewed Plate . 


For the sake . of completeness (since the flat plate is 
the limiting case of the helicoidal shell being considered) 
and due to the relatively small amount of attention given 
the skewed cantilever plate in the literature, plate eigen- 
values and mode shapes for various aspect ratios and skew 
angles will be presented. Comparisons of the present solu- 
tion with other analytical and experimental results will be 
given for the skewed plate. 

Four aspect ratios (AR = z/2b ) are considered: 0.5, 1.0, 

2.0, 4.0. Poisson’s ratio is set to a nominal value of 0.3 
for all plate calculations. Also a is restricted to 0°, 15°, 

30°, and 45°. Permissible ranges are 0 < AR < ~ and 0° <_ a < 90°. 
The authors have made unverified calculations for a up to 85° 
without any difficulty. However, those results will not be 
included since the helicoidal shell data was restricted to 
0° _< a <_ 45°. The values chosen for AR and a yield 1 6 possible 
parameter combinations. All calculations used a 35 term double 
power series for w. No loss of numerical accuracy (16 place 
arithmetic was used) was apparent for the 35 x 35 matrix eigen- 
value problem. As a increases, the coupling between odd and 
even n terms increases, causing the computation of the eigen- 
values to take longer. The total cost of the 16 eigenvalue 
problems was $20.00 on the IBM 370 digital computer. Conver- 
gence curves will not be given since only the 5 lowest modes 
of each parameter combination will be given. For 35 terms, 
the 5 lowest modes are converged (throughout the text the word 
"converged" means that less than 0.5$ change occurred in the 
eigenvalues when additional terms were taken in the double 
power series). 

Table 1 contains the 5 lowest eigenvalues for all 1 6 
parameter combinations. For the data given in Table 1, eigen- 
value curves may be plotted versus a, since the eigenvalues 
are continuous functions of a. Plots of the eigenvalues versus 
a for AR = 0.5, 1.0, 2.0, and 4.0 appear in Figs. 5, 6, 7, and 
8 respectively. For AR = 0.5 and 4.0 (Figs. 5 and 8) no mode 
crossing can be detected for 0° <_ a <_ 45° although nodal 
patterns are altered with increasing a (see Figs. 9, 10, 11, 
and 12). For AR = 1.0, Fig. 6 shows the 3rd and 4th mode 
shapes reversing their order at very nearly a = 32°. Also 
the 5th and 6th mode shapes reverse position at a = 38°. Thus 
one frequency must generate 2 independent mode shapes at points 
of intersection of 2 curves (i.e., at some precise a value). 

These reversals of position in the frequency spectrum of a 
given mode shape can be evidenced in the w nodal patterns 
(Figs. 9-12). For AR = 2.0 only the 5th and 6th mode shapes 
cross (nearly at a = 41°). As aspect ratio increases, a given 
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nodal pattern may drastically change its position in the 
frequency spectrum (see Pigs, 9-12), This is due to the 
fact that, as AR decreases the relative longitudinal stiff- 
ness increases, thus causing more of the lower modes to have 
longitudinal nodal lines. 

Tables 2 and 3 contain comparisons of these power 
series results to experimental and other analytical results. 
Aspect ratio is 1.0 in both tables. 

a - 30° (Table 2) 

Column 1 contains experimental results from References 
[11] and [12]. Column 2 contains the authors' direct solu- 
tion. Column 3 has analytical results due to Claassen [13], 
who used 18 beam functions in a Ritz analysis. Note that 
Claassen' s results are all higher than the author’s direct 
solution. Since the Ritz method yields upper bounds, the 
author's numbers are more converged than Claassen' s [13] re- 
sults. Column 4 . contains experimental results given in 
Reference [14]. The agreement between analysis and experi- 
ment is very good considering that in some instances, shakers 
were used to excite the plate (shaker position can drasti- 
cally affect mode shapes and thus nodal patterns, as evidenced 
in Ref. [14]). 

a = 45° (Table 3) 

Columns 1, 2, and 4 are the same as in Table 2. Column 
3 contains analytical results obtained using Reissner's 
theorem [15]. In light of the experimental data in Columns 
1 and 4, and the well converged (less than 0.5# error) analy- 
tical results in Column 2, some of the results in Column 3 
are upper bounds and some are lower bounds. Again the 
analytical results in Column 2 agree favorably with the experi- 
mental results. 
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THE HELICOIDAL SHELL 


Since there is no data for comparison to the cantilever 
helicoidal shell, convergence of the eigenvalues must be 
checked to establish their validity. In choosing the three 
double power series, no preference is made for u, v, or w. 

All three double power series are given the same number and 
degree terms. For large aspect ratio higher order £ terms 
are needed, whereas for smaller aspect ratios higher order n 
terms are needed. 1 6 place arithmetic is used throughout 
the calculations. For the largest number of terms taken (30 
terms in, each series; a 90 x 90 matrix eigenvalue problem) 
no loss of accuracy was noted. 

First a simple configuration (&/R = 0.2, h/R = 0.01, 
a = 0°, and y - 0.1 rad) was checked. For the helicoidal 
shell aspect ratio is AR = £/2Ry. In this particular case 
AR = 1.0. Convergence curves for the first 5 eigenvalues 
appear in Fig. 13. The first 4 modes, are converged after 
only 16 terms (48 x 48 matrix) and all 5 are converged at 24 
terms (72 x 72 matrix). Next, the same shell with a = 45° 
was checked. The convergence curves are given in Fig. 14. 

For 24 terms all 5 eigenvalues are very close (within 1.5%) 
to the converged numbers of 30 terms. As a last example, an 
aspect ratio 4.0 shell (l/R = 0.8, h/R = 0.01, a = 15° and 
y = 0.1 rad) was checked. Convergence curves are shown in 
Fig. 15. Again 24 terms yield results within 2% of the con- 
verged numbers of 30 terms. Convergence of other examples, 
such as small aspect ratio (large y ), have been, checked but 
will not be included, since the results are virtually the 
same as the three cases presented. Of course it would be 
desirable to run all cases with 30 terms to assure converged 
numbers. For a 24 term run (72 x 72 matrix), the average 
cost on the IBM 370 digital computer was $10.00. Thus, as 
a matter of economy, and in order to present the widest range 
of data possible, all results presented will be for 24 or 25 
terms. As has been shown, the eigenvalues for 24 terms are 
within 2% of the converged numbers of 30 terms. 

The geometric parameters are 2./R, h/R, a and y. Poisson's 
ratio was set to a nominal value of 0.3 for all calculations. 
For thin shell theory h/R is usually taken in the range 
0<h/R<0.05. This is in keeping with the approximations made 
in the derivation of the thin shell equations. In this paper 
only h/R in the range 0 . 005<h/R_<0 . 05 will be considered. 

Also a will be restricted to 0_<a£ir/4 and y to 0.1 rad <y<0.7 
rad. Permissible ranges for a and y are of course 0<_a<w72 
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and 0 < y £ ir. The permissible range for a/B. is 0 < A./R < °°. 
But again, in order to limit the enormous amount of possible 
parameter combinations, l/R is taken in the range 0 , ±<_l/R<_ 0 . 8 . 
In order to observe, the effects of a given parameter, all 
other parameters are fixed at values which could represent a 
typical shell. 


Effect of h/R 
For a - 0°'; &/R = 0.2; y = '0.1 rad 

The aspect ratio for this case is 1,0, For small R/h 
(= 20), the primarily bending modes (denoted by B]_, B 2 ,,,,) 
have relatively high frequencies. Also one encounters prima- 
rily membrane modes (denoted by Mi, M2,...) much earlier in 
the spectrum. But as R/h increases, the frequencies of the 
primarily bending modes decreases rapidly (see Fig. 16). 

Since the frequencies of the primarily membrane modes change 
very little with R/h, they appear later and later in the 
spectrum due to the rapid decrease in the bending mode eigen- 
values. 


For ct = 30 °, &/R = 0 . 2 ; y = 0.1 rad 

The trend of the eigenvalues versus R/h remains the same 
with this change in a (see Fig. 17). The eigenvalue for the 
Mi mode is consistantly lower for a = 30° than for a = 0°. 
Also the eigenvalues for B4 and B 5 are consistantly lower for 
a = 30° than for a = 0°. For a = 30°, modes B3 and B4 are 
very close in frequency versus R/h, although the mode shapes 
differ considerably. 


Effect of a 

For l/R = 0.2; h/R = 0.01; y = 0.1 rad 

The aspect ratio is 1.0. Table 4 contains the first 5 
eigenvalues for a = 0°, 15°, 30°, and 45°. Also included are 
the relative amounts of bending strain energy and membrane 
strain energy, and the generalized mass. The eigenfunctions 
were normalized in such a way as to make the generalized mass 
be between 1.0 and 10.0. The generalized stiffness (total 
strain energy) is then Xg X generalized mass. The generalized 
mass and stiffness would be used in calculating transient 
motion of the shell by the normal coordinate technique. If a 
mode has more than 60 % of the strain energy due to membrane 
energy, it is considered a primarily membrane mode, denoted by 
M. If more than 60% of the strain energy is due to bending, 
it is considered a primarily bending mode, denoted by B. 
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If neither is more than 60 % , the mode is denoted by B-M or 
M-B, the first letter indicating which is greater in magni- 
tude. Eigenvalues will not be graphed versus the remaining 
geometric parameters. This is because mode shapes and fre- 
quencies change rapidly and. it becomes very difficult to draw 
curves for the 5 lowest modes. Rather, the 5 lowest eigen- 
values will be given in tabular form along with figures 
showing the w nodal patterns. For the present case, w nodal 
patterns appear in Fig. 16. The 1st and 2nd mode shapes can 
be easily followed versus a. The 3rd and 4th modes, however, 
cross between a = 30° and 45°. Also between 30° and 45°, a 
bending mode crosses a membrane mode (5th mode). 


For ft/R = 0.8; h/R = 0.01; y = 0.1 rad 

The aspect ratio is 4.0. Eigenvalues and energy percent- 
ages are given in Table 5. In Fig. 19 the first 3 modes can 
easily be followed versus a. However, the 4th and 5th modes 
are difficult to follow. At a = 15°, the 4th and 5th mode w 
nodal patterns are very nearly the same. Yet, the relative 
magnitudes and/or shapes of u and v may be entirely differ- 
ent. This fact will be illustrated in a later section. After 
a = 15° , the 4th and 5th modes seem to retain their position 
in the frequency spectrum. 


Effect of y 

For &/R = 0.8, h/R = 0.01; y = 0° 

The aspect ratio varies with y (AR = £/2Ry). Eigenvalues 
and strain energy percentages are given in Table 6. w nodal 
patterns appear in Fig. 20. For y = 0.1 rad, the w nodal 
patterns are very similar to those for the flat plate. As y 
increases, the relative longitudinal stiffness increases caus- 
ing two things to occur: first, more longitudinal node lines 
appear in the lower modes, and second, the membrane energy for 
•the lower modes increases considerably. It is almost impossible 
to follow a given nodal pattern versus y* Although one or two 
deflection components may have the same (or similar) nodal 
patterns for different nodes, the third component may have a 
very different nodal pattern for the two modes. This is 
illustrated by modes 2 and 3 for y = 0.7 rad-. Even though the 
w nodal patterns are different, u and v have similar nodal 
patterns (see Fig. 21). The relative magnitudes of u, v, and 
w are also similar for the two modes shown in Fig. 21. Upon 
close inspection, although u and v nodal lines are similar, 
their positioning is different. 
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For l/R = 0.8; h/R = 0.01; a =30° 

Aspect ratio varies with y as noted before. Eigenvalues 
and strain energy percentages are given in Table 7« w nodal 
patterns are shown in Fig. 22. Modes 1, 2, and 3 seem to 
follow an orderly progression from y = 0,1, 0.3» 0.5, to 0.7, 
with each increase in y bringing an additional longitudinal node 
line. As previously mentioned, this is a result of the rela- 
tive increase of the longitudinal stiffness. Modes 4 and 5 are 
difficult to follow. Mode 4 seems to follow the trend of in- 
creasing number of longitudinal node lines with increasing y. 

At y = 0.5 rad, modes. 4 and 5 have similar w nodal patterns; 
but, u and v nodal patterns may be quite different. 


Effect of l/R 
For h/R = 0.01; a = 0; y = 0.1 rad 

The aspect ratio (AR = l/2Ry) varies since l/R varies. 
Eigenvalues and strain energy percentages are shown in Table 
10. w nodal patterns are shown in Fig. 23. Unfortunately, 
Fig. 23 is somewhat distorted. y is fixed at 0.1 rad, so as 
l/R decreases, the length of the planforms should decrease. 
Instead, the length of the planform has been fixed at one 
unit and the width (2yR/0 thus increases as l/R decreases. 
Except for one membrane mode each in l/R = 0.8 and 0.4, all 
other modes shown are primarily bending modes. As AR de- 
creases ( l/R decreasing), the relative longitudinal stiff- 
ness increases, with the result that more of the lower modes 
have longitudinal node lines. Notice the similarity in 
Fig. 23 and the analogous flat plate in Fig, 9. The w nodal 
patterns are identical with the exceptions of mode 5 for 
AR = 2.0 and modes 4 and 5 for AR = 4.0. 


For h/R - 0.01; a = 30°; y = 0.1 rad 

Again the aspect ratio varies with l/R. Eigenvalues and 
strain energy percentages are given in Table 9- w nodal 
patterns are shown in Fig. 24. Fig. 24 has the same distort- 
ions as Fig. 23. The 1st and 2nd modes can be followed easily 
as AR decreases ( l/R decreasing). As before, more of the lower 
modes have longitudinal nodal lines as AR decreases. Note the 
similarity to the flat plate nodal patterns (Fig. 11). The 
reason for the similarity is the shallowness (y = 0.1 rad) of 
this particular shell. 
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CONCLUDING REMARKS 


Free vibration modes and frequencies have been presented 
for both the cantilever skewed plate and cantilever helicoidal 
shell with variations of all geometric parameters. From the 
convergence curves in the case of the shell, and from experi- 
mental and other analytical results in the case of the plate, 
the authors feel that the results are numerically sound. 

The shell results point up the serious problems which 
confront the designer. That is the question of what the first 
few mode shapes are. It is impossible (especially for large y 
in the case of the shell) to guess at the mode shapes in order 
to get a 1 or 2 term Rayleigh-Ritz solution that will yield 
anything close to the true fundamental frequency. In short, 
anticipating the mode shapes of the helicoidal shell is not 
feasible. As R/h increases, the spacing between shell eigen- 
values rapidly decreases, causing mode excitation, identifica- 
tion, and resolution to be very difficult for the experimenta- 
list. Also as a and y increase, nodal patterns change drasti- 
cally, thus furthering the problem of mode identification. At 
some a values two mode shapes exist for one frequency, for both 
the plate and shell. Positions of given nodal patterns in the 
frequency spectrum are altered drastically with changes in 
aspect ratio for the plate and shell. 

It is felt that this paper represents a third step toward 
the analytical solution to the flutter of the thermally stressed 
turbine blade. The first and second steps are represented in 
Refs. 2 and 16 respectively. The turbine blade, immersed in 
a non-steady fluid flow field and a non- steady structural 
acceleration field, poses a challenging problem for an analy- 
tical solution. The Ritz method, when power series are pro- 
perly employed with finite time intervals and finite space 
intervals in Hamilton’s Law of Varying Action, appears to offer 
a direct method of attack on this and other most difficult 
problems . 
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APPENDIX I 


COMPARISON TO EXACT SOLUTIONS FOR SPHERICAL CAPS 
AND THE COMPLETE CONICAL SHELL 


A literature search was made including the Applied 
Mechanics Reviews , Index to Contemporary Acoustical Literature 
published by the Acoustical Society of America, a computer 
survey from NASA Langley Research Center, and the NASA mono- 
graphs, "Vibration of Plates" and "Vibration of Shells" by 
A. W. Leissa Cl], [173. No exact solution exists for the 
skewed cantilever plate or for the cantilevered helicoidal 
shell [173. No numerical results from either calculation or 
experiment, were found in the literature on the vibration 
modes and frequencies of the cantilevered helicoidal shell. 
However, comparison of calculations by the authors, using power 
series, have been made to the "exact" solutions for two other 
shell geometries and are given below. 


The Deep Spherical Shell 

The exact solution for the deep spherical shell is given 
in terms of Legendre functions of the first kind by Kalnins 
[18] and Kraus [19]. In the exact solution, (h/R)2/i2 has 
been neglected in comparison to unity, in order to arrive at 
a solvable differential equation (this is in keeping with the 
approximations made in thin shell theory). In the present 
analysis (h/R)2/l2 was retained. 

For axisymmetric motion (i.e., only u and w deflections 
present), Kalnins [18] has given eigenvalues for various 
boundary conditions. In the direct solution by the authors, 
a 10 term power series was used for both u and w. A compari- 
son of results for a clamped base and a free base is given in 
Table 1 of this Appendix. The agreement is quite evident. 

The reason the direct solution (shown by Ritz [M] to be an 
upper bound for stationary problems) yields a lower eigenvalue 
for the first few modes than the exact solution, is felt to be 
that the differential equations that would be obtained from 
the energy equation are slightly different than the different- 
ial equations of Ref. [17] due to the neglect of (h/R)2/i2 in 
comparison to unity. 
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TABLE 1. DEEP SPHERICAL SHELL 



Axisymmetric Motion 


h / R = 0.05 , v = 0.3 , = 60° 

X* = u Rs/p/E ' 

Mode No. 

Free Base 

Clamped Base 

1 

o.93i a v> 

1.006 

[0.9308]° 

C1.0056] 


1.088 

1.391 

2 

[1.088] 

[1.391] 


1.533 

1.697 


[1.5336] 

[1.6963] 

A 

2.348 

2.375 

*+ 

02.3488] 

[2.3746] 


2.544 

3.486 

O 

[2.5478] 

[3.4872] 

Pi 

3.497 

3.991 

o 

[3.5346] 

[3.9907] 

7 

4.951 

4.974 

[5.0956] 

[4.9792] 


a. unbracketed numbers due to Kalnins [18] 

b. [ ^bracketed numbers from author’s direct 


solution 
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The Complete Cone 


Due to the geometric singularities at the apex, the com- 
plete conical shell ha!s received little treatment by the Ritz 
method, and limited exact solution. Neglecting tangential 
inertia (u = 0), an exact solution has been obtained by using 
Airy’s stress function and the method of Probenius [20]. 

Dreher [20] used the Donne 1 1- Mu shtari shell equations. The 
author used Novozhilov’s strain-displacement relations to 
obtain the direct solution of the free base case, with and 
without tangential inertia* Reference [20] gave numerical 
results for the clamped and free base, A comparison of results 
is given in Table 2 of this Appendix. The small difference in 
_ the frequencies from the exact solution and the direct solu- 
tion, when tangential inertia was neglected, is probably due 
to the fact that slightly different equations were solved. 

The result of neglecting tangential inertia has the predictable 
effect of raising the natural frequencies since mass is effec- 
tively taken out of the system (or conversely, constraints 
have been added) . 


29 


TABLE 2. CONICAL SHELL 



Axisymmetric Motion 
24.25 p = 0.000254 Ib-sec/in 4 

v = °'3 0 = 60° 

E=10 r psi N # = 12(1 -Z'*)lf/’h)/tan z /3 

h = 0 • 4626 in =1 0 4 

Frequency , Hz 

N A rvHo 

Neglecting u 

Including u 

No. 

Exact Solution, 
Ref. [20] 

Direct Solution, 
author 

Direct Solution, 
author 

1 

912.36 

912.35 

891.47 

2 

1322.53 

1322.60 

1269.34 

3 

1823.71 

1828.55 

1746.54 

4 

2414.92 

2410.31 

2310.01 

5 

3122.80 

3264.53 

2987.97 
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APPENDIX II 


MATRIX ELEMENTS 


n ijkjt = o 


t&i 


for Ll±jQ 
odd 


even 


I2 ljki ■ 0 


o 1 .li 
TO 1 +R 


odd 


even 


I3i -fif t = 2 (k+l) J+l(l+l) 
( l+k+2 ) (j+Jt ) 


odd 


l4 iJkA * 2 — * 4#l .. 

( I+k+2)( J+i) 


even 


odd 


even 


I 5 i Jkl - 2 


odd 


even 


16 


ijki 


o 1 .11 

'( i'+k+3 ) 7?F T 


odd 


even 
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for (.1+1) 

I7 ijkl * 0 odd 

( J+4+1) even 

l8 ijki - 2 fi+t+vy 13^17 odd 

m 0 even 

I9 uki * 2 jr+l+n) f$+f^T odd 

* 0 even 

«°uki - 2 -^n 0<1<J 

“ 0 even 

I11 ijkl " 0 odd 

" 2 (i+ScH') ( J+/-TJ * Ten 

I12 ii« - 2 Ws) 1? T3+TT odd 

■ 0 even 

H3ljicl * 0 odd 

- 2 T^Hrf^T 
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114 


ljk| 


g U+k+ ST IM+ H 


for U±l) 
odd 


even 


I15 1Jk , - 0 


- g 


odd 


even 


Il6 lJkjL " 0 


= 2 


(MfMr 


odd 


even 


i i7 uki - ° 


odd 




even 


IlOlJki - 0 


= 2 




odd 

even 


rai»» ‘ g 

= 0 even 

I2° ljk , - 2 ( Ii| + j 1+I j (1*8) Ham 

odd 

= 0 even 
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Mlijkl * 0 


for (J+A) odd 


" 2 (t+w3T(M+vi for U+J) eTCn 

M2 ijki * M1 ijki 


M3ijkl * Il4 iJkJi 


K1 iJkA “ ^ijkt + (i/Rr) 2 (sin 2 a + 

-(sin a)(A/RY)i 3 lt j kl 

KSiJia - (*/*v)( il§ 2 Llfcos a)l 4 1Jkl + v(co» a)I 5 tJkl 
-(l/RyXsln a)(cos a) ( 1 | v ) I 6 ij k j) 

K 3 ijki * (i/Ry)v(cos a)( 17 ^^ - CA/Ry)(sin a) • 
l8 ijki ) 

K4 ijkl - {(l+(h/R) 2 / 12 )(A(cos a)/R 7 ) 2 + (l+(h/R) 2 / 3 ) 
• IlgSl (I(Bin a)/R 7 ) 2 } I 2 1Jkl 

+ (l+(h/R) 2 / 3 )li^lI Il 1JkA - (A/R<y) • 

(sin a) I3ijki 
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K5ijki ■ (l(cos a) A) 2 *8±3*il/y “ ((h/R) 2 /(X2y)) * 
{(lAY) 2 (l+(l-v)sin 2 a) I9ijki + v I10 ijkl 

- 2(iAr)(sin a) Illijki + 2(l-v) I12 1Jkl 

- 2(l-v)(i/RY)( sln a ) I^ijk*} 

K6 1Jki - (*/* ) 2 H4 1JklL (eos 2 a) + ((h/R) 2 /12) • 

{(R/i(cos a)) 2 I15ijkl + (iAY 2 (cos a)) 2 ' 
H6ij k ji + (tan 2 a + v) I17ijk£/Y 2 
+ (2(l-v)+4tan 2 a) IlSijj^/y 2 - 2 (^/Ry) • 
(tan a) I19i jkl/Y 2 ( c ° 8 a) - 2(tan a) • 
(R/£y) I20 1;Jkl /(cos a)} 

KBijkt “ I1 5ijkl+ (*/b) 4 Ilfitjkl + (sin 2 a+vcos 2 a)» 
(i/b) 2 I17ijki + (2(l-v)cos 2 cw-4sin 2 a) • 
(i/b) 2 Il8 1Jk£ - 2(i/b)3(sin a) I19ijkJL 

- 2(i/b) (sin a) I20 1Jk ^ 
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PLATE EIGENVALUES 



83.7205 


78.3589 


67.5967 































TABLE 2. SKEWED 


a =30°; v-. 

Mode 


No. 

Experiment , 
Refs. |11], [ 12 J 

Direct Solution, 
Be res 

1 

2.87 

2.9522 

2 

6.92 

7.0878 

3 

18.38 

19. 0494 

4 

19.15 

19.4758 


• 



5 


30.48 


31. 1485 
























TABLE 3. SKEWED PLATE EIGENVALUES 


a =45°; v = 0-3 ; AR =1.0 

Mode 

No. 

X p = wVph/D ‘ /cos 2 a 

Experiment , 
Refs. [11], [12] 

Direct Solution, 
Beres 

Relssner's Theorem, 
Ref. [15] 

Experiment, 
Ref. [i 4 ] 

1 

2.13 

2. 2688 

2.06 

2.31 

2 

5.53 

5.6992 

5.63 

5.48 

3 

13.26 

13.6363 

13.86 

13.31 

4 

15.06 

16.0858 


15.90 

5 

25.09 

25.7958 



























TABLE 4. HELICOI 


J/RzO-2 h/RxO-OI 

a =0° 


Mode 

No. 

Vp(1- itVE/COScC 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

1 

0.055033 

16.99 

83.01 

9.1242 

2 

0. 12315 

0.4O 

99.60 

1.6942 

3 

0.31920 

6.17 

93.83 

9- 0570 

4 

0. 39306 

0.79 

99.21 

6.5230 

5 

0.45040 

0.63 

99-37 

3.0463 


a =30° 


Mode 

No. 

X5 - Cj * 

VpO-vVE' /cosec 

°lo of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

1 

0.051386 

7.86 

92.14 

5.4171 

2 

0. 12108 

4.62 

95.38 

6.7129 

3 

0.32675 

2.29 

97.71 

4.6053 

4 

0.33108 

1.63 

98.37 

4 . 5698 

5 

0.53356 

10.22 

89.78 

1.5491 


I 


SHELL EIGENVALUES 


v =0-3 T'O-I rad 



a =45° 






TABLE 5. HELICOIDAL SHELL EIGENVALUES 
i/R =0-8 h/Rs001 v =0-3 r= 0-1 rad 


Mode X S ' <x>* 
No. 



% of strain % of strain generalized 

energy due energy due mass 

to membrane to bending 


% of strain % of strain generalized 

energy due energy due mass 

to membrane to bending 



No. ■ 



°/o of strain 
energy due 
to bending 

generalized 

mass 

9?.27 

1.9527 

90.58 

2.8079 1 

89.52 

6.6308 

54.00 

2.4239 

53.05 

8.2090 


Mode X s - co* ° fo *°f strain % of strain generalized 
No energy due energy due mass 

Vp(1- A'EVcosec to membrane to bending 


1 0.009890 


2 0.050025 


3 0.087251 


0. 15069 


0. 20298 







































































































TABLE 6. 


HELICOIDAL SHELL EIGENVALUES 


Jt/ R = 0-8 h/R = 001 v =0-3 a=0° 

7=0-3 rad 


7 =0- 5 rad 


Mode 

NO. 

A S = «• 

Vp(1- v'VE'/coscc 

°/o of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

a 

0.049641 

69.01 

30.99 

2.6410 

2 

0.061378 

18.17 

81.83 

1.3951 

3 

0,14914 

•89.96 

10.04 

1.5511 

4 

0. 17898 

33.19 

66.81 

9.4340 

5 

0.17959 

39.35 

60.65 

9.5614 


Mode 

NO. 

A s = co ■ 

Vp(1- v^)/E7cOSoC 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

n 

0.053481 

48.87 

51.13 

2.1032 

2 

0.066911 

56.84 

43.16 

6.3?77 

3 

0. 13920 

53.49 

46.51 


a 

0. 18983 

36.29 

63.71 

2.8389 

5 

0.20016 

61.72 

38.28 

1.6372 


Mode 

No. 

A s = co • 

p /o of strain 
energy due 
to membrane 

°lo of strain 
energy due 
to bending 

generalized 

mass 

Vp(1- I'VE/cosec 

n 


18.11 

81.89 

2.1815 

2 

0.052883 

89.74 

10.26 

1.3872 

3 

0. 15929 

30.59 

69.41 

3.6941 

a 

0. 15974 

8.51 

91.49 

1.5465 

5 

0.20723 

38.07 

61.93 

3.7834 


7=0-7 rad 


7=0-1 rad 


Mode 

No. 

Vpfl- 7)/E7 cosoc 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

m 

0.013817 

20.41 

79.59 

ns— 

2 

0.085751 

19.30 

80.70 

2.1607 

3 

0.099644 

0.33 

99-67 


D 

0.23240 

99.91 

0.09 

5.9594 

5 

0.24289 

16.69 

83.31 

2.6617 

























































































TABLE 7. HELICOIDAl SHELL EIGENVALUES 


A/R = 0-8 h/R = 001 „=0-3 a =30° 

7=0-1 rad 7=0-3 rad 


Mode 

No. 

A 5 ” CO • 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

Vp( 1 - vVE'/coscc 

i 

0.019676 

21.59 

70.41 

1.0371 

j 2 

0.069371 

50.44 

49.56 

1.7562 

3 

0 . 12*106 

31.37 

68.63 

1.7545 

D 

0.15527 

28.02 

71.98 

3.8605 

5 

0.24472 

31.03 

68.97 

3.8958 


Mode 

No. 

As - CO • 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

Vp(1- K*)/E'/cosoc 

■ 

0.011659 

7.73 

92.27 

1.9527 

2 

0.064690 

9.42 

90.58 

2.8079 

3 

0.10389 

10.48 

89.52 

6.6308 

a 

0.18166 

46.00 

54.00 

2.4239 

5 

0.22456 

46.95 

53.05 

0.2090 


7=0- 5 rad 7=0-7 rad 


Mode 

No. 

A s = CO- 
Vp(1-^)/E7cosoc 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

H 

0.025434 

32.31 

67.69 

1.1078 

2 

0.094740 

47.90 

52.02 

7.1253 

3 

0.13482 

68.49 

31.51 

1.6844 

a 

0.14516 

44.22 

55.78 

1.6861 

5 

0.20398 

66.03 

33.17 

1.6626 


Mode 

No. 

A s = co • 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

Vp(1- v/)/E7cosoc 

1 

0.023129 

23.54 

76.46 

1.6050 

2 

0.082180 

42.44 

57.56 

3.5571 

3 1 

0.13052 

32.89 

1 

67.11 

6.0257 

a 

0.14811 

1 

83.05 

16.95 

2.5619 

5 

0.16823 

39.59 

60.41 

2.2914 

































































































TABLE 8. HELICOIDAL SHELL EIGENVALUES 


h/R = 0-01 v =0-3 a=0° 7=0-1 rad 


X/ R.0-1 z//R = 0-2 


Mode 

No. 

Vp(Vv^)/E7co5oc 

°/o of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

1 

0.055033 

16.99 

83.01 

9* 12*12 

2 

0.12315 

0.40 

99.60 

1.6942 

3 

0.31920 

6.17 

93.83 

9.0570 

4 

0. 39806 

0.79 

99.21 

6.5230 

5 

0.45040 

0.63 

99.37 

3.0463 


Mode 

NO. 

A s = CJ- 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

Vp(1- v^/E'/cosoc 

1 

0. 10619 

9.75 

90.25 

1.1958 

2 

0.15516 

0.67 

99.33 

1.0221 

3 

0. 29847 

0.49 

99.51 

8.8034 

4 

0.56742 

0.15 | 

99-85 

3.5975 

5 

0.63881 

1.19 

98.81 

1.6759 


^/R = 04 ^/R = 08 


Mode 

No. 

Vp(1- 7)/E7cosec 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

1 

0.02769 

19.53 

80.47 

2. 1162 

2 

0.10712 

0.29 

99.71 

1.6425 

3 

0. 16707 

15.00 

84.92 

2.3493 

4 

0.35020 

0.49 

99.51 

2.7004 

5 

0.41555 

99.97 

0.03 

2.2949 


Mode 

No. 

X s = <-)■ 

% of strain 
energy due 
to membrane 

% of .strain 
energy due 
to bending 

generalized 

mass 

Vp(1- v/)/E7cosoc 

1 

0.013817 

20.41 

79.59 

6.4158 

2 

0.085751 

19.30 

80.70 

2. 1607 

3 

0. 099644 

0.33 

99.67 

3-9252 

4 

0. 23240 

99.91 

0.09 

5.9594 

5 

0.24209 

16.69 

83.31 

2.6617 





TABLE 9. HELICOIDAL SHELL EIGENVALUES 


! h/R = 001 V =0-3 a =30° 7=0-1 rad 

J/RrnO-1 j(/R-0-2 


Mode 

No. 

Aj = CO • 
Vp(1- v»Ve7cosoc 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

D 

0.051386 

7.86 

92.14 

5.4171 

B 

0.12108 

4.62 

95.38 

6.7129 

B 

0.32675 

2.29 

97.71 

4.6053 

B 

0.33108 

1.63 

98.37 

4.5698 

B 

0.53356 

10.22 

89.78 

1.5491 


Mode 

NO. 

As = to ■ 
Vp(1- v'VE'/coscc 

°/o of strain 
energy due 
to membrane 

■°lo of strain 
« energy due 
to bending 

generalized 

mass 

B 

0. 10377 

3. 80 

96.20 

1.3865 

B 

0.15990 

4.03 

95.97 

9.5907 

3 

0. 28535 

0.40 

99.60 

2.4998 

B 

0.52231 

0.16 

99-84 

1.5343 

5 

0.65908 

i 1 

1.07 

98.93 

2.0466 


Jt/R = 0-4 .7/R =0- 8 


Mode 

NO. 

A 5 - to ■ 
VpCl- vVE'icosoc 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 


0.011659 

7.73 

92.27 

1.9527 


0.064690 

9.42 

90.58 

2.8079 

3 

0.10389 

10.48 

89.52 

6.6308 

4 

0.18166 

46.00 

54.00 

2.4239 

5 

0.22456 

46.95 

53.05 

8.2090 


Mode 

No.. 

A 5 - CO • 

% of strain 
energy due 
to membrane 

% of strain 
energy due 
to bending 

generalized 

mass 

Ypd- */)/E '/cosec 

B 

0. 024325 

7.86 

92.14 

. 1.5070 

B 

0. 10082 

5.09 

94.91 

4.6684 

3 

0.16408 

7.83 

92.17 

1.5731 

B 

0.29471 

1 

9.75 

90.25 

3.8707 

5 

0.35614 

92.90 

7.80 

3.2429 
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FIGURE 4. RIGHT CIRCULAR 
CYLINDER CUT ALONG y = ir 
AND DEVELOPED INTO THE 
PLANE , SHOWING THE HELI - 
COIDAL SHELL AND SKEWED 
COORDINATE SYSTEM 
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c , degrees 

SKEWED PLATE 
[S vs a, AR=0.5 
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FIGURE 6. SKEWED PLATE 
EIGENVALUES vs cc,AR=1.0 






FIGURE 8. SKEWED PLATE 
EIGENVALUES vs cc,AR=4-0 
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Figure 9. SKEWED PLATE , w NODAL PATTERNS 
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FIGURE 13. CONVERGENCE OF EIGENVALUES 
FOR HELICOIDAL SHELL 
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'ERGENCE OF EIGENVALUES 
.ICOIDAL SHELL 
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i/R=0-8 , h/R = 0-01 , a=15°, 7 = 0-1 rad 
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FIGURE 15. CONVERGENCE OF EIGENVALUES 
FOR HELICOIDAL SHELL 










































FIGURE 21 . DISPLACEMENT NODAL PATTERNS FOR 
HELICOIDAL SHELL, i/R = 0-8, h/R=0 01, a=0°, 7=0-7 






















